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SUMMARY 


An independent variable different from the time for elliptic orbit integration 
is used here. Such a time transformation provides an analytical step-size regu 
lation along the orbit. An intermediate anomaly (an anomaly intermediate be- 
tween the eccentric and the true anomaly) is suggested for optimum perform- 
ances. A particular case of an intermediate anomaly (the elliptic anomaly) is 
defined, and its relation with the other anomalies is developed. 


INTRODUCTION 

The development of integration schemes offering an automatic step-size control 
(such as the Runge-Kutta-Fehlberg scheme (ref. 1) or the variable-step vari- 
able-order multistep method developed by Shampine and Gordon (ief. 2)) may be 
an invitation to return to the direct integration of the differential equations 
of the orbital motion in Cartesian coordinates with time as the independent var- 
iable. If the automatic step-size control does take care of a step-size reduc- 
tion near pericenter and a step-size enlargement near apocenter, the step-size 
distribution along the orbit is not as smooth as through an analytical step- 
size regulation. This is because the step size is predicted on the basis of a 
current local truncation error estimation. So that this step size is not re- 
jected too often during the following step, a conservative factor is intro- 
duced, prohibiting an excessive predicted step size. As a result, the integra- 
tion accuracy along the orbit is not homogeneous. 

Analytical step-size regulation means that a new independent variable t is 
chosen (often called fictitious time) so that equal steps in x produce the de- 
sirable distribution of time steps along the orbit. The time is no longer the 
independent variable, but it becomes a new dependent variable to be integrated 
together with the other state variables. The drawback of having one more equa- 
tion to integrate is compensated for by a more efficient use of the function 
evaluations. In addition, a highly efficient fixed-step, multistep integration 
method can be used, which will be free of the overhead involved in a variable- 
step method. 


THE GENERALIZED SUNDMAN TRANSFORMATION 

To achieve an analytical step-size regulation for the integration of orbits, a 
change in the independent variable 

dt = fdx (2.1) 

is proposed, where t is the physical time, x the new independent variable 
and f a function of the state. Equation (2.1) is called a time transforma- 
tion. By choosing f in the form 

f = (2.2) 
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where r is the radius vector, a a real number and a coefficient, the 

new independent variable can be of an angular type describing 2 tt over one rev- 
olution like the classical anomalies used in celestial mechanics. Equation 
(2.2) is referred to as the generalized Sundman transformation. 

Case 1 : a = 0 

By choosing 


Cq =Vla|3/p 

a is the semimajor axis, ]x is the universal gravitational constant times the 
mass of the central body, and t is equal to the mean anomaly. 

Case 2 ; a = 1 

Assuming an elliptic orbit, by differentiating Kepler's equation 


t = - e sin u) 


(2.3) 


where t is the time since pericenter pass and u the eccentric anomaly, and 
by using the relation 

r = a(1 - e cos u) . (2.4) 

equation (2.1) combined with equation (2.2) becomes 

y-1/2 C(x“^ _ g ggg u)^”°^ du = dt (2.5) 


where the time is replaced by the eccentric anomaly. With a = i , equation 
(2.5) becomes 

p-1/2 ci“^ a^/2 du = dT 


By choosing 

c-| =V"a/i7 

du = dT and T is the eccentric anomaly. 
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Therefore, equation (2.2) becomes 


f 



r 


which is the original Sundman transformation. 

Case 3 : a = 2 

The term function of u in equation (2.5) can be expressed in terms of the 
true anomaly 4> as 


du 

■y1-e^(l - e cos u) 


d<l> 


Therefore, by choosing 


1 


C2 “ 

yyad - e^) 

r2 



defines a time transformation where the new independent variable is the true 
anomaly. 


THE INTERMEDIATE ANOMALY 

At this stage, it is interesting to examine the distribution of steps along an 
elliptic orbit defined by the three transformations of the last section. 

Figure 1(a) shows how 12 equal steps in mean anomaly are distributed along an 
0.8 eccentricity ellipse. All steps are accumulated around the apocenter, and 
the pericenter is completely depleted. This is the result of using an indepen- 
dent variable proportional to the time. 

Figure 1(b) shows how 12 equal steps in eccentric anomaly are distributed. A 
significant improvement over the preceding case is noticed, but steps around 
pericenter are still sparse. 

Figure 1(c) shows a distribution of 12 equal steps with respect to the true 
anomaly. This time, a sparse number of steps around apocenter is noticed. 
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By simple inspection of figures 1(b) and 1(c), it can be anticipated that a 
more adequate distribution of steps along the ellipse would be reached by 
having a time transformation defining an independent variable whose characteris 
tic is intermediate between the eccentric and the true anomaly. 

The concept of intermediate anomaly is introduced here by a pure qualitative ar 
gument. It will be shown now that one particular intermediate anomaly is the 
mathematically logical extension to the ellipse of the polar angle for the cir- 
cle. 


THE ELLIPTIC ANOMALY 

Recalling that the eccentric anomaly corresponds to Of = 1 of equation (2.2) 
and the true anomaly corresponds to d = 2 , an intermediate anomaly is defined 
for any value of ct such that 1 < a < 2. As suggested by Nacozy (ref. 3) the 
first obvious choice is a = 1.5. This case will be investigated in more de- 
tail as follows: 


Case 4; a = 3/2 


Equation (2.5) becomes 


du 


C 3/2 yT - e cos u 


= dx 


(4.1) 


If the auxiliary angle 6 = ^/2 is introduced, equation (4.1) becomes 


O.5C3/2 + e) = 


d6 


yi - sin^ ® 


(4.2) 


with 


k 2 = 2 e /(1 + e) 

1 

The right-hand side of equation (4.2) is precisely the argument of the normal 
elliptic integral of the first kind 



Therefore, equation (^.2) becomes 


2F(e, k) 

T = - 

C 3/2 u( 1 + e) 


(4.3) 


where T is measured from pericenter. Equation (4.3) ^ direct relation 

ship between the true anomaly <j) = 26 and the intermediate anomaly t. As an 
elliptic function is involved in its definition, this particular choice of an 
intermediate anomaly will naturally be called elliptic anomaly and will be 
denoted from now on by v. 

At apocenter, cj) = it and 


F(iT/2, k) = K(k) 


where 


K(k) 



sin^ 9 


is the complete elliptic integral of the first kind. 

The constant G 3/2 is chosen such that v = it at apocenter, therefore 

2K(k) 

C 3/2 = — ■■ ' 

TT'yyd + e) 

and equation (4.3) becomes simply 


TrF(6, k) 

V = 

K(k) 


(4.4) 


Equation (2.2) is now defined as 


2K(k) r3/2 

IT -^pd + e) 


(4.5) 
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Figure 1(d) shows the distribution of 12 steps equally spaced with respect to 
the elliptic anomaly. In contrast to the cases of the eccentric and true 
anomaly, the distribution of steps given by the elliptic anomaly favors neither 
the pericenter nor the apocenter and seems to be ideally suited for an 
analytical step regulation. 


THE ELLIPTIC ANOMALY AS ARGUMENT OF THE 
JACOBIAN ELLIPTIC FUNCTIONS 

In order to express the elliptic anomaly in terms of the traditional anomalies, 
it is necessary to introduce the inverse of the elliptical integral of the 
first kind: the Jacobian elliptic functions (ref. 4) 


sn(v, k) = sin 6 

cn(v, k) = cos 6 (5.1) 


such that 


V = F(6, k) 


The functions sn and cn are periodic of period 4K(k) and are a generaliza 
tion of the circular functions sine and cosine for the ellipse. For zero 
eccentricity (k = 0) 


sn(v, 0) = sin v 
cn(v, 0) = cos V 

By introducing equation (4.4) into equation (5.1) 

sin 0 = sn(K V/1T, k) 
cos 6 = cn(K v/tt, k) 

the true anomaly 4» = 6/2 is expressed in terms of the elliptic anomaly, as 
follows: 


sin (}) = 2sn(n, k) cn(n, k) 
cos 4) = 1 - 2sn^(n, k) 
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where 


n = K v/iT 


By using the usual expressions relating the true anomaly to the eccentric 
anomaly, the eccentric anomaly u can be expressed directly in terms of the el 
liptic anomaly as follows: 


sin u 



sn(n, k) cn(n, k) 
1 - k^sn^Cn, k) 


cos u 


1 + e - 2sn.2(n, k) 

1 + e - Besn^Cn, k) 


Equation (2.4) expressing the radius r becomes (in terms of the elliptic 
anomaly) 


a(1 - e) 

r = ■■ ■ 

1 - k2sn2(n, k) 


By using the focal unit vectors i directed toward pericenter and normal 
to i., position r and velocity v are expressed as 


r = r(1 - 2sn2(n, k))i + 2r sn(n, k) cn(n, k)j 


( r - i)2sn(n, k)cn(n, k) 

, a(1 - e2) 

+ (1 - 2sn^(n, k))J dn(n, k) 

Formulas involving the elliptic anomaly were first developed in reference 3, 
where an intermediate variable proportional to the elliptic anomaly is defined. 

The three time transformations corresponding to the eccentric, true, and ellip- 
tic anomalies were already outlined in reference 5 and are developed along an 
alternate way in reference 6. 


V = yjvJ 1+e)/r 
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THE NUMERICAL EVALUATION OF ELLIPTIC FUNCTIONS 


An apparent drawback to using the elliptic anomaly as an independent variable 
is the need to estimate elliptical functions at each integration step. 

However, fast and accurate algorithms for computing elliptic integrals and func- 
tions are given in reference 7 • 

The algorithm for estimating the Jacobian elliptic functions sn and cn is 
based on the Gauss transformation and iteratively applied. A total of only two 
trigonometric functions is required. For an accuracy of at least 15 decimal 
places, no more than 100 elementary operations are needed. 

If just the time transformation equation (eq. (4.5)) is used, only the complete^^ 
elliptic integral of the first kind K(k) is required, and the corresponding 
algorithm does not involve more than 20 elementary operations. 


CONCLUDING REMARKS 

The use of an independent variable different from the time leading to an 
analytical step-size regulation is the most efficient way to compute elliptic 
orbits. 

Instead of the commonly used eccentric or true anomaly as a new independent var- 
iable, it is suggested to use an intermediate anomaly (intermediate between the 
eccentric and the true anomaly). One choice among the possible intermediate 
anomalies is the elliptic anomaly, so called because it involves elliptic 
functions. - The elliptic anomaly appears to be the most natural way for 
extending the polar angle of the circle to the ellipse. 

The anomalies discussed in this paper were used as the independent variable dur- 
ing a numerical comparison. An orbit of eccentricity, 0.73, was chosen, and a 
time element was included. Results of the comparison showed that the use of 
the elliptic anomaly as the independent variable led to an accuracy about one 
order of magnitude higher than any other classical anomaly. These results will 
be presented in a subsequent paper. 
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(c) The true anomaly 
a = 2 



(d) The elliptic anomaly 
a = 3/2 


Figure 1.- Concluded. 
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